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Spacecraft  Charging  Model- 
Two  Maxwellian  Approximation 


I.  INTRODUCTION 

In  Garrett  *  and  Garrett  and  DeForest,  “  it  was  demonstrated  that  the  plasma 
at  geosynchronous  orbit  is  to  first  order  representable  by  expansions  in  terms  of 
Maxwellian  distributions.  In  particular,  a  description  in  terms  of  the  sum  of  two 
Maxwellian  components  was  found  to  adequately  represent  the  plasma  in  about 

HO  percent  of  the  cases  studied.  In  Garrett,  3  a  simple  spacecraft  charging  code 

*  4  5 

based  on  the  work  of  Whipple  and  DeForest'  was  developed  which  made  use  of 

the  discrete  particle  measurements  returned  by  the  ATS- 5  geosynchronous  satel¬ 
lite  in  the  energy  range  50  eV  to  50  keV.  The  model  suffered,  however,  from 
being  unable  to  treat  accurately  Maxwellian  components  except  in  a  narrow  range 
of  potentials  and  temperatures  due  to  the  necessity  oT  taking  discrete  energy  steps 
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m  performing  integrations  (that  is,  a  temperature  of  10  eV  would  require  at  least 
1  eV  resolution;  a  temperature  of  lit)  keV  would  then  require  over  40,  000  such 
steps  if  both  were  to  he  resolved  simultaneously).  It  became  clear  that  the  only 
solution  t « >  such  a  problem  was  to  carry  out  the  detailed  calculation  of  the  currents 
associated  with  a  Maxwellian  distribution  function. 

In  the  first  part  of  this  re[uirt  we  w  ill  review  the  formulation  developed  in 
Garrett.'^  As  each  term  is  discussed,  the  necessary  integral  will  be  developed 
and  the  results  for  a  single  Maxwellian  presented  for  various  temperatures  and 
potentials.  The  final  section  will  review  the  details  of  the  two  Maxwellian  model 
and  conclude  with  a  comparison  between  theory  and  actual  observations.  These 
results  indicate  a  standard  deviation  of  tl,  300  V  between  observation  and  theory 
over  the  range  -300  V  to  -10,000  V,  while  the  region  of  applicability  of  the  model 
is  between  *100  V  and  -30,000  V  for  potentials  and  between  10  eV  and  a0,  000  keV 
for  plasma  temperatures  —  much  larger  than  that  of  Garrett. 


2.  CKMKII.  rm-om 


The  fundamental  problem  in  determining  the  potential  on  a  surface  immersed 
in  a  plasma  is  the  balance  of  currents  leaving  and  striking  the  surface  —  an  im¬ 
balance  always  results  in  a  net  addition  or  subtraction  of  charge,  altering  the 
potential  until  equilibrium  is  reached.  In  this  study  we  will  consider  only  the  fol¬ 
lowing  currents  lor  a  single  point  (that  is,  leakage  currents,  capacitance,  etc., 
are  ignored): 


J, 


(J, 


JSE+ JSI  *  JBSE  +  JPh'  ° 


(1) 


where 


JE 


incident  electron  current 
incident  ion  current 


JSE 


secondary  emitted  electron  current  due  to  Jp 


JSI 


secondary  emitted  electron  current  due  to  .1 


USE 


backscattered  electron  current  due  to  Jp 


pn 


photoelectron  emission  (independent  of  Jp  or  J  j> 


Given  the  incident  ion  and  electron  spectra  at  the  satellite  surface,  the  cur¬ 
rents  Jp,  Jj,  Jgp,  J^j,  and  (Jpn  will  be  ignored  in  this  study,  though  it  is 

implicitly  included  in  the  model  and  can  be  readily  modeled;  see  in  particular 


Garrett  und  IV1  orest  >  are  round  and  the  potential  p  on  the  spacecraft  varied 
until  Kq.  1 1)  holds. 

Subsequent  sections  will  outline  the  procedure  for  <  aleulating  each  of  the  cur¬ 
rents  in  Kq.  (1)  as  functions  f  N.  the  number  densitv,  T,  tin’  Maxwellian  tem¬ 
perature,  and  o,  the  satellite  t,  snare  potential.  W  e  will  m  fact  be  dealing  with 
•be  nunibei  flux.  Nl  .in  units  of  parti  le  >  rm”-sei  -st  rather  than  the  actual 
current  density.  1  or  a  truly  omnidirectional  flux  they  are  related  by: 

J  q  *  (NFV  (2) 

where  q  electronic  charge. 

Alsi  ac  will  be  making  the  limiting  assumption  of  a  thick  sheath  surrounding 
lb<  satelliti  .  n.at  is,  w e  will  assume  that  the  distribution  function  ?  at  the  sur- 
'■i>  e  of  the  space  raft  is  related  to  the  ambient  distribution  function  f.  for  a 
Maxwellian  distribution,  by 

?<K't  e'q,>  K1  f( F: •  >  (3) 

w  here 

.•  •>  _  1/  U  T 

ftn  mm  jtk n  -  c 

and 

K  Boltzmann’s  constant 

M  mass 

1  energy  in  ambient  medium 

K'  energy  at  surface 

K  K'  *  q  •  o  (4) 

At  geosvtu  ht  onous  orbit,  the  sheath  (that  is.  the  region  over  which  particle 
orbits  arc  distorted)  is  much  larger  than  the  spacecraft  (typical  values  are  <0  m 
for  the  sheath  as  compared  w  ith  a  ft  w  meters  for  the  satellite)  so  that  Kq.  (3)  is 
a  good  approximation.  At  lower  iltitudos  (particularly  in  the  plasmasphere)  this 
relation  breaks  down  and  the  detailed  geometry  of  the  satellite  and  its  sheath  must 
he  considered. 

f’.  Barrett,  II.  B.  and  Def  orest.  S.  K.  (11'71'a)  Time-varying  photoeleet ron  flux 
effects  on  spacecraft  potential  at  geosynchronous  orbit.  .1.  Geophvs.  lies. 
84:283-21  -  - 
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In  terms  of  the  incident  particle  distribution,  f,  at  the  surface  of  the  satellite, 
we  can  write  the  equations  for  the  incident  electron  and  ion  currents  as  follows: 


J  i:ft i.i  di: 
0 


(5) 


Substituting  I'q.  (3)  t  r  !',  and  rrmpinbering  that  the  ambient  particles  do  not  have 
energies  below  0,  we  find  for  the  electrons: 


-N,.  /  Mj,  V  e  -a:-0>  t 

«>,)  /  ’ 


For  0  ■>  0,  x  0.  implying: 


h:  kknftk"(1  '*  ri:' 


loro  -  0.  x  0,  implying: 


KlAi: 1  !  ' 


(0  te) 

e 


iv  here 

Mj.  mass  of  electron 

3 

N  number  density  o f  electrons  in  n  cm 
Tj,  electron  temperature  in  cV 
0  satellite  potential  (in  units  of  eV) 

K  energy  of  electrons  in  eV 

Kj.  >.  <2(1  X  1 0 ' '  for  Jj.  in  n  cm- -sec-sr. 


(6a) 


(lib) 


(6c) 


Throughout,  we  are  using  the  number  flux  in  the  equations  even  though  we  refer 
to  it  as  the  current  —  the  two  are  related  by  llq.  (2).  Also  the  units  of  charge  q 
will  be  incorporated  in  0,  so  that  0  is  in  the  units  of  eV  rather  than  V  (the  s i gn 
of  the  charge  will  be  explicit  in  each  equation).  rf.  and  Tj  are  assumed  to 
include  K.  the  Boltzmann  constant. 
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electron  number  FLUX 

8- 


LEGEM 

<P  >  ♦  iOv 
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Fi^re  la.  Number  Flux  of  Electrons  (Jed  in  Units  of  n/cm2-sec-sr  as  a  Func 
tn>n  of  Maxwell -Boltzmann  Temperature  Tt^)  for  Various  Values  of  Satellite 
Potential  (0).  N ,,  is  assumed  to  be  1  cm "2 


ION  NUMBER  Flux 

9~ 

8- 


LEGEND 

0  =  +  10V 
0  >  0V 
0  -  100V 
0  •  -  IOOOV 
0  ■  -  5000V 
0  *  -  10,000V 


8  10  12 
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Eigutc  lb.  Number  I  lux  of  Ions  (Jr)  in  Units  of  n  cm'-sec-sr  as  a  Function  of 
Maxwell-Boltzmann  Temperature  (T j>  for  Various  Values  of  Satellite  Potential 
(0).  Nj  is  assumed  to  be  1  cm'-’ 
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Similarly,  for  the  ions  (assumed  to  be  protons): 


For  p  <  0,  x  -<J,  implying: 

Jj  KjNjTj  2  (l-o  Tj) 

For  ij>  ?■  0,  x  0,  implying: 


(7a) 


(7b) 


where 

Nj  number  density  of  ions  in  n/cm' 

T.  ion  temperature  in  eV 

5  2 

Kj  1.2437X10  for  Jj  in  n/cm  -sec-sr. 

Expressions  (fi)  and  (7)  are  plotted  in  Figure  1.  As  we  shall  see,  they  are  of 
fundamental  importance  in  many  spacecraft  charging  applications  and  can  be 
easily  applied  to  calculate  the  current  to  the  satellite  due  to  the  ambient  environ¬ 
ment  for  a  plasma  consisting  of  two  or  more  >  irnponents.  As  an  example,  for  the 
case  of  a  two-Maxwellian  distribution  of  electrons  and  a  negative  potential: 


J,,  Kj,.  (N1E 


TlEl/2e*/T1E  +  N2E 


T2E1/2  e*/T2E) 


where 

N1E  number  density  for  electron  component  1 
TIE  temperature  for  electron  component  1 
N2E  number  density  for  electron  component  2 
T2E  temperature  for  electron  component  2 


(8) 


Similar  equations  follow  for  combinations  of  3,  4,  or  higher  sums  of  Maxwellian 
components. 


1.  SECONDARY  EMISSION  Cl  RKK.NT 


Secondary  emission  is  the  result  of  electrons  being  elected  from  the  satellite 
surface  due  to  impinging  electrons  and  ions  (ion  secondaries  are  not  important). 
The  secondary  emission  current  is  given  by  (remember,  and  J^  are  in  units 

of  n/cm2 -sec-sr): 
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.1  I  f.  (1  'U 1  > I  .11 

M'  *■' 


...  d-\  t  !\t  s.  .  mnlai  v  olei  t  roil  viol. I  is  appr. 'Xlmnteil  for  I  lie  ele.-t  i  .mis  hv : 


K<U>«  FI'''1  •  e-"1’*'  n 


lor  0  *  0 


loll-  1  IV 


for  o  *  o 


.»  10.1' 


(  l .  }l  l  i 

O 

n.ii  ’ •  l  !Mio.  i 

I  In  lim.  non  ,v  ii  I,  for  0  •  0,  has  Ih'imi  pintle. 1  m  I'ijjitro  11  was  rhosett  to  be 
in  .  I.  ■  .i  >;  i  .  i  i  .nt  \v  it  li  t  In-  sii'otuln  r\  loot  i  on  vie  hi  for  oloot  ions  ini  part  nifi  on 
aliiioiiiuin  prosenl ml  in  Whipple. 

I  or  oloolrotlS  .mil  C  ■  f>|.  U'l  l'i  oolites. 


.1  /  .>.,.(1  '  f..(Kir  .11' 

f-i  J  I 


K  •  l  \  or,.  /  _(K  1  °K  •*'  ,  -1 

‘  n  •  ,  •  e  '  /  •'  1  <>■•’*  .')•  I'  -IK 

1  I'  1 


I'ho  a  bo  vo  integral  oannot  bo  intO|iratoil  nnatvt  ieall\ .  11  o  liavo,  however,  n 
I I'^jriit oil  it  mini i' ri i'ii  1 1 \  vising  the  linuss-l  aiiiiorro  ipiiulraturo  fornnilns.  I  lio 
results,  tint-mall  .ed  to  .1  T  .,  if>).  are  plotte.l  in  l  ipiii  r  :l.  1'  is  dear  from  tins 

figure  that,  for  tin  nui^r  shown  (nanielv,  10  oV  ■  I',.  ■  100  ke\  an. I 

It  .  0  -  10  keV): 


cndary  electron  yield 


(12) 


JSE<TE'*^  S(Te) 

where 


r 

-0.  028(1TV 

<t>  >  0 

S(Tr.) 

\  °- 

(1411 

-0.842c  fc‘ 

T.,  <  200  eV 

p  «■  0 

\ 

-0.  OOOUT., 

L*. 

V.  0. 

i4:i 

+  0.  740  e 

T,,  ?  200  eV 

$  £  0 

Although  we  have  not  used  this  expression  in  the  model  presented  here,  the 
gain  in  computation  time  from  this  latter  expression  should  outweigh  any  loss  in 
accuracy.  As,  however,  we  have  sought  accuracy  in  this  study,  we  have  in  every 
case  integrated  Eq.  (11)  as  indicated. 

The  expression  for  the  secondary  emission  of  electrons  due  to  ions  was 
taken  to  be: 

'  K2  e-U/<E+V> 

6 j(t:>  (  Kj 

(° 

w  he  re 

K.  0.  oat;  ,  K.,  5 

U  •10(10  ,  V  200  ,  K  700 

This  is  plotted  in  Figure  4  over  a  plot  of  6,(E)  for  aluminum  taken  from 

■1  1 

Whipple.  As  the  deviation  above  100  keV  is  significant,  we  have  limited  our 
model  to  50  keV  (well  above  the  region  of  naturally  occurring  charging  phenomena), 
so  the  deviation  is  not  of  concern  in  our  calculations. 

Inserting  Eq.  (El)  into  Kq.  (0)  we  find,  for  <j)  5  0: 


<  0  E>ll 

^  <  0  H  >  E  (ED 

d  >  0 


JSI  f  6,(E)TI(E)E  dE 

si  M“  J  l  T,  ~ 


K.N.  ~d>  /  T . 


1  -d> 


R 


f  ”E/Ti  f 

S  J  e  E  dE  +  K2  J 


K.  j  e 


c'K/Tl.  e-U/(E+V)EdE 


(14) 


where  Mj  is  the  mass  of  the  ions. 
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5.  IlM'.kSC.  Vn'K.RKl)  Kl.K  TRON  UKKINT 


As  implied  in  its  name,  the  backscattered  current  is  the  result  of  the 
backseattering  or  reflection  of  the  incident  particles.  II  is,  for  our  calculations, 
only  a  significant  factor  for  electrons.  The  basic  equation  is: 


BSE 


OC 

—  /  dE’  J  B(E\  E)Fj,,(E)E  dE 


( 171 


ME  x 


where 


lO  if  <t>  *  0 

\o  >f  <f>  *  0 

H(K',  El  percentage  of  electrons  scattered  at  a  given  energy  K'  as  a 
result  of  an  incident  electron  at  energy  E. 

E'  energy  of  backscattered  particles  (E1  *  1C) 

E  energy  of  incident  particles 

Follow  ing  Sternglass,  ‘  IKE',  E)  is  usually  given  in  the  normalized  form: 


nil''  K) 

IKE',  E) 


(18) 


G  is  given  as  a  function  of  E'/E  for  aluminum  by: 


G(E'/E)  ~  *  (E'/EKl  -  E'/E) 


(10) 


Wr  obtained  Eq.  (19)  from  an  approximate  fit  to  Sternglass  and  De Forest 
(private  communication,  see  Garrett'S.  Substituting  into  Eq.  (17): 

Fop  <&  ^  0: 


JBSE 


(20a) 


7.  Sternglass,  E.J.  (1954)  Barkscattering  of  kilovolt  electrons  from  solids, 
Phys.  lies.  95:345. 


18 


For  o  ^  0: 


USE 


11.  <1  «  »  T  )  t-  1.2  •  -  L3 


(3  +  2<I,/Te) 


-v 

_ - 

(y  +  <?  tf) 


<i> 


(20b) 


where: 

1 

nete 

L1  f 

nete 

N 

T3/2 

1  E 

•> 

N 

\J  p 

-i.3  /2 

The  integration  is  performed  by  the  Gauss -l.aguerre  quadrature  formula. 
The  results,  normalized  by  Jj,(d,  Tg),  are  plotted  in  Figure  0.  The  detailed 
numerical  integral  results  in  Figure  6  can  be  roughly  fitted  by: 


JRSE  BS<*,  Te)  •  JE<*.  Te) 


(21) 


where 


j  <t>  s0 

BS(d,  T.J  /  n 

'  (0.253  -  0.0027  e0,0'^)  +■  (0.0807  +  0.657  e'0’  184°)  . 


-(0.00017  t  0.0008  p-0-231^  te 


0^0 


As  before,  the  actual  integral  was  used  for  the  calculation. 
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1  igure  < ■ .  Hatio  of  Jlaokscatterod  Electron  Nunl't  I  Flux 
Due  to  Electron  llombardment  t.l  nsi  ’  1,1  Ambient  Electron 
Number  Flux  at  Surface  (Jj.O  Versus  Temperature  ( )' j.O 
for  Various  Potentials 


i 


(>.  r»t>  MWMI  l  l\\  MUtiVl 


As  discussed  in  the  Introduction,  Garrett1  and  Garrett  and  Do  Forest*1  have 
presented  evidence  that,  at  least  for  the  energy  range  -  10  eV  to  '  100  keV,  the 
plasma  population  can  be  adequately  represented  b\  .i  sum  of  two  Maxwellian 
components  for  both  the  electrons  and  ions.  Experience  shows  this  t,  be  true  for 
approximately  80  percent  of  the  cases.  Although  more  components  would  naturall 
give  a  better  fit,  and  in  about  10  percent  of  the  cases  a  single  Maxwellian  is  ade¬ 
quate.  a  two  Maxwellian  approximation  is  a  very  convenient  way  to  represent  the 
plasma  as  it  can  be  directly  derived  from  the  first  four  moments  of  the  distribu¬ 
tion  function.  These  four  moments  (see  Garrett)  are  the  quantities  most  often 
provided  by  particle  experiments,  and  hence  most  readily  available.  With  these 
factors  in  mind,  our  charging  model  was  designed  to  be  optimized  for  two 
Maxwellian  dist ribut ions. 


Kq.  <11  becomes: 


IJK(0.  Nli:.TtK)  *  J  |.(o,  N2E,  T2K1I  -  ||  J  ,(o,  N  11.  T 1 0  *  .1  to,  N21,  T21>| 

'  Ust;<0.  mi:,  T1K)  ♦  JSK(o.  N2K,  T2K)|  *  [J S1(C.  Nil,  I'  1 1)  +  Js](o,  N21,  T21) I 
*  lJBSK(d*  NlE*  T1K'  ‘  ‘'hsK10-  X~1;-  T2K,I  f  0  <22> 


w here: 


Ml',  I'll:  Maxwellian  number  density  and  temperature  for  electron 
component  1 

N-’K,  T2K  Maxwellian  number  density  and  temperature  for  electron 
component  2 


N2I,  121 


JK 


Nil,  I'll  Maxwellian  number  density  and  temperature  for  ion 
component  1 

Maxwellian  number  density  and  temperature  for  ton 
component  2 

incident  electron  current  given  by  Kq.  (it) 

incident  ion  current  given  by  Kq.  (71 

secondary  electron  current  given  by  Kq.  (Ill 

secondary  ion  current  given  by  Kq.  (la) 

backscattered  electron  current  given  by  Kq.  (211 

photoelectron  current,  assumed  zero  in  this  paper  (see  Garrett,' 
however) 


S  K 
JSl 


BSK 

JPH 


Kq.  (221  is  solved  by  substituting  values  of  o  until  it  is  satisfied,  a  process 
for  which  a  variety  of  schemes  can  be  employed.  In  this  study,  o  was  steadily 
decreased  in  coarse  steps  until  Kq.  (221  became  negative.  .At  this  point,  the 
interval  bracketed  by  tins  value  and  the  preceding  value  was  subdivided  into  finer 
intervals  and  the  process  iterated  until  a  close  approximation  to  zero  was  found. 

He  fore  analyzing  the  results  of  our  model,  one  final  point  remains  to  be  dis¬ 
cussed  concerning  the  adequacy  of  our  assumptions.  The  assumption  of  a  thick 
sheath,  or,  equivalently,  assuming  the  satellite  to  be  a  point,  limits  the  accuracy 
of  our  results  as  the  sheath  can  alter  the  ambient  plasma  to  a  significant  degree 
before  the  plasma  strikes  the  satellite.  This  assumption,  however,  is  not  serious 
in  comparison  with  that  of  assuming  the  satellite  to  be  made  uniformly  of  aluminum. 
First,  the  material  properties  (that  is,  6  and  Ill  are  at  best  only  qualitatively 
understood.  Secondly,  it  is  obvious  that  few  satellites  consist  of  aluminum. 
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In  fact,  the  twi>  geosynchronous  satellites  ATS-  >  and  ATS-t>  considered  in 
this  study  consisted  of  a  variety  of  materials.  In  order  to  compensate  for  this, 
in  Garrett'*  a  careful  analysis  was  carried  out  in  order  to  study  the  effects  of 
varvittg  6  and  E.  It  was  found  adequate  to  multiply  these  factors  by  a  constant  f 
the  order  of  unity  to  bring  the  results  for  ATS-5  and  ATS-G  into  agreement  with 
observations.  The  correction  factors  adopted  for  tins  study,  based  on  ATS-5, 
w  e  re : 

6,,  1.3  by 

^  0.  55  6j  (2d) 

ft(E',  E)  0.  4  B(E',  E) 

The  ATS-t>  values  were  different,  but  as  we  are  interested  in  a  "universal"  model 
wc  will  employ  the  ATS-  >  values  throughout  as  representative. 

T.  1*  \T  \  WU.VSIS 

ATS-5  and  ATS-fi  are  both  geosynchronous  satellites.  We  have  integrated  the 
digital  data  from  the  Vniversitv  of  California  plasma  detectors  on  these  two  satel¬ 
lites  to  obtain  the  four  plasma  moments  (for  the  energy  ranges  50  eV  to  50  keV 
and  0  to  80  keV.  respectively)  for  the  ambient  environment.  To  be  precise,  the 
ambient  environment  immediately  preceding  eclipse  entry  (or  following  eclipse 
exit)  and  immediately  after  exit  (or  immediately  preceding  exit)  from  eclipse  was 
computed,  the  results  being  corrected  for  satellite  potential.  The  corresponding 
eclipse  potentials  were  also  computed.  The  iwo  Maxwellian  components  where 
possible  (if  not,  a  single  Maxwellian  was  chosen)  were  computed.  These  values 
and  their  averages  (giving  three  sets  of  2d  values  each)  were. employing  the  model 
as  rust  outlined,  used  to  compute  the  potentials  for'  the  sunlit,  eclipsed,  and  aver¬ 
age  test  cases.  The  three  sets  of  estimated  potentials  are  plotted  versus  the 
observed  potentials  in  Figures  7a,  7b,  and  7c.  Although  there  is  obvious  agree¬ 
ment  between  the  observed  and  predicted  potentials  in  these  figures,  we  will  demon¬ 
strate  mathematically  their  agreement  in  the  subsequent  paragraphs  (see  Table  1). 


Table  1.  Statistical  Analysis  ot  A  I  S-  >  and  A  I  S-li  I  .  lipse  Data 


Number  of  Point s  23 
Degrees  of  freedom  21 
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The  first  set  is  for  sunlit  data,  eclipse  data  make  up  the  second  set,  and  tin 
third  set  is  an  average  of  the  two. 
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OBSERVED  POTENTIAL  UV> 

a.  Sunlit  Data 

Figure  7.  Observed  Versus  Predicted  Potentials  for  A'I'S-.i 
and  ATS-ti  for  Three  Sets  of  Pstimates  of  the  Ambient  Pnvi- 
ronmetit  —  Sunlit  Conditions,  Kclipse  Conditions,  and  the 
Average  of  Tin  sc  Two  Pstimates 


c.  Average  Data 

Figure  7.  Observed  Versus  Predicted  Potentials  for  ATS-'> 
and  ATS-fi  for  Three  Sets  of  Estimates  of  the  Ambient  Envi¬ 
ronment  —  Sunlit  Conditions,  Eclipse  Conditions,  and  the 
Average  of  These  Two  Estimates  (Continued) 
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a  'i  \iMir \i  \n\i,>si> 


It  there  is  a  relationship  between  the  experimental  and  theoretical  values  of 
the  potential,  the  simplest  model  should  be  a  straight  line  relation  of  the  form 

V  A  •  BX  (24) 

w  here  X  signifies  the  experimental  values  and  Y  the  theoretical  values  of  the 
potential.  For  perfect  agreement  between  X  and  5  we  should  have  A  0  and  B  1. 

I  hree  sets  of  data,  each  having  23  points,  have  been  obtained  and  each  set 
has  been  studied  statistically.  For  the  first  set,  the  least  squares  estimates  b  and 
a  of  B  and  A,  respectively,  are  found  to  be 

b  0.714  and  a  202.  (iO  . 

lo  create  a  .  onfidenee  interval  for  A,  we  choose  o  0.05  as  the  level  of  signifi¬ 
cance,  and  we  conclude  that  a  !>f>  percent  confidence  interval  for  A  is 

Ml  <  A  < M2 

with  Ml  a  -  tp  Q25  an;i  M2  a  •  t^  jy,.  being  the  lower  and  upper  con¬ 
fidence  limits,  respectively.  For  the  data,  these  were  found  to  be  Ml  -654.98 
and  M2  10’>0.  US.  I  his  shows  that  9  a  percent  of  the  time  the  mean  A  of  the 
5-  intercept  of  the  straight  lint'  (24)  will  fall  inside  the  above  confidence  interval. 

In  other  words,  the  confidence  interval  would  contain  the  parameter  A  of  the  popu¬ 
lation  with  a  probability  of  95  percent. 

As  mentioned  before,  t.  have  a  perfect  relationship  between  X  and  Y  we 
should  halt’  A  0  and  R  1.  Now,  the  question  is  asked,  is  it  logical  to  consider 
the  value  of  A  to  be  zero  or  not”  For  this  we  test  the  null  hypothesis  HO:  A  0, 
against  the  alternative  HI:  A  \  0.  Here  we  consider  a  two-tailed  critical  region 
for  a  0.0a  level  of  significance.  The  region  of  acceptance  would  be  all  values  of 
a  such  that 
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^  w \ 1 1  accept  the  hypothesis  tf  the  estimated  value  a  ol  A  is  inside  the  region 
of  acceptance  and  reject  it  if  it  is  in  the  critical  region.  For  A  0  tin-  region  of 
acceptance  becomes 

-tl.r»7.  58  <  a  «  857.  58  . 


Since  a  2011.80,  it  falls  inside  the  acceptance  region  and  the  null  hypothesis 
(A  0)  is  accepted.  This  procedure  guarantees  that  the  probability  of  rejecting 
the  hypothesis  A  0  when  it  is  in  fact  true  is  less  than  »  0.05.  In  addition, 
this  procedure  guarantees  the  probability  of  accepting  A  0  when  it  is  not  true  is 

smaller  than  in  any  other  test  with  the  same  it.  (This  is  the  so-called  uniformly 
most  powerful  test  when  the  alternative  hypothesis  is  A  '  0). 

A  similar  situation  prevails  for  the  slope  B  of  the  line  (24).  A  (1  -  <>)100% 
confidence  interval  for  the  parameter  It  is  given  by  the  double  inequality 


h  *  SI3  *,>  2 


b  i  S  t 

U  t> 


which  means  that 

,,|i  [h  -Sn‘„/2  '  B  '  b  '  Sb‘d/2]  1 

t  boosing,  as  before,  n  0.  05  as  the  level  of  significance',  we  obtain  the  following 
confidence  interval  for  H: 


N 1  <  H  «  N2 


where  N1  0.4!*  and  N2  0.04.  This  shows  that  05  percent  of  the  time  the  slope  11 
of  the  straight  line  (24)  would  he  between  N  1  and  N2. 

Now  testing  the  hypothesis  II0:  11  1  against  the  alternative  HI;  It  *  l  we 

consider,  as  before,  a  two-tailed  critical  region  for  ■>  0.05  level  of  significance. 

The  region  of  acceptance  would  be  all  values  of  h  such  that 


I’ll 


0.  Of 


or  all  values  of  b  which  satisfy  t In'  double  inequality 


K1  <  b  ■.  K2 


where  K1  11  -  S.^  I 


and  Is 2  ]l  l  S.|  t 
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In  thi'  present  case.  K1  0.777,  K2  1.223,  and  h  0.714.  It  appeals  that  in 
tins  ease  tin'  null  hypothesis  H  1  is  not  consistent  with  the  data,  and  thus  our 
host  estimate  (minimum  variance)  would  he  H  0.714.  However,  it  should  he 
noted  that  for  a  level  of  significance  0.01  considering  a  critical  region  outside 
of  the  interval 


0.097  s  It  ^  i.;t02$) 
the  hypothesis  It  l  is  accepted. 

I'ahle  1  gives  a  general  picture  of  the  statistical  analysis  done  for  all  three 
sets  of  data.  The  letter  A  signifies  the  V-intereept  of  the  line  (24).  H  denotes  the 
slope,  and  SA  and  S((  represent  the  standard  deviations  of  A  and  It,  respectively. 
IIP  denotes  the  null  hypothesis.  111  the  alternative  hypothesis.  Die  standard  error 
is  denoted  hy  o  and  the  correlation  coefficient  is  denoted  by  r.  Hased  on  this  table 
it  appears  that  to  a  reasonable  degree  of  accuracy  (tl,  200  V),  our  model  based  on 
a  tvvo-Maxivellian  approximation  and  a  simple  material  configuration  predicts 
spnnvraft  potentials  in  eclipse. 


*>.  COM  I  I  SION 


111  the  preceding  sections  we  have  outlined  a  procedure  for  taking  a  two- 
Maxwellian  distribution  (or,  for  that  matter,  any  number  of  Maxwellian  distribu¬ 
tions)  and  computing  the  potential  difference  between  a  shadowed  satellite  surface 
and  the  ambient  environment.  Our  results  indicate  an  accuracy  (actually  the 
standard  error  between  the  observed  and  predicted  potentials)  of  tl,  200  V.  Phis 
is  to  be  compared  with  the  "exact"  calculation  using  the  detailed  distribution 
function  of  Oarrett'  for  which  an  error  of  t-700  V  was  reported  and  that  of 
f  .nrrett  et  al  tor  the  crude  estimate  of  if  =  -Tj,  I  ,n  (.1^,  10J  )  of  t2,  000  V.  As  the 
range  ol  accuracy  of  the  representation  of  the  distribution  function  and  current 
balance  decreases  from  Garrett'  to  this  paper  to  Garrett  et  al,  ,‘1  this  trend  is  in 
good  agreement  with  expectations. 

1  lie  power  ol  the  techniques  presented  in  this  paper,  however,  is  not  the 
accuracy  (which  can  be  improved  by  "tuning"  the  model  as  discussed  in  llarrett'S 
but  in  the  possibility  of  deriving  an  efficient,  analytic  expression  for  the  potential. 
Hus  is  evident  in  Kqs.  (ti),  (7).  (12).  (1(1).  and  (21).  Multiple  Maxwellian 
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distributions  ran  bo  oasil\  inoorporatod  using  thoso  rquat ions  in  a  straightforward 
fashion.  Wo  intond  to  pursuo  tins  formulation  m  subsoquont  papors,  loading 
ultunntolv  to  :»  simplo  mot  hod  of  doming  spaooo  ra  ft  lc  s|m»v  potentials. 
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